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Abstract 

In 1942, the famous astronomer and physicist Belgian Georges Lemaitre (1894-1966) published 
in the Bulletin of the class of Sciences of the Royal Academy of Belgium, an article in which he 
proposes a method that he calls "rational iteration". He rediscovers both the Aitken process 
and the Steffensen’s method. 

We will see how the scientihc concerns of Georges Lemaitre in the year 1942 (resolutions of 
differential equations) have led to this rediscovery. 

Key words Lemaitre, Aitken, Steffensen, A^, hxed point, ordinary differential equation, 
Picard iterative process. 


Resume 

En 1942, le celebre astronome et physicien beige Georges Lemaitre (1894-1966) publie 
dans le Bulletin de la classe des sciences de I’Academie Royale de Belgique, un article 
[Lemaitre G., 1942 a dans lequel il formule une methode qu’il nomme iteration rationnelle. 
II redecouvre a la fois le processus du A^ de Aitken et la methode de Steffensen. G’est 
en 1927 que le mathematicien neo-zelandais Alexander Graig Aitken (1895-1967) publie son 
fameux procede d’acceleration de convergence [Aitken A. G., 1927) . Quelques annees plus tard, 
en 1933, le mathematicien danois Johan Frederik Steffensen (1873-1961) expose sa methode 
[Steffensen, J. F., 1933[ . La methode de Steffensen est vue d’ailleurs a prfeent comme un 
corollaire du A^ de Aitken0. 

Si le compte-rendu de Particle de Georges Lemaitre paru dans le Jahrbuch de 1942 mentionne 
la methode de Steffensen[§, Lemaitre ne fait pas reference aux deux mathematiciens. Gela pent 


*Institut de Mathematiques de Bourgogne, Universite de Bourgogne, France, leferran@u-bourgogne.fr 

1. Steffensen ne mentionne pas Aitken dans son article. Le compte-rendu paru dans Zentralblatt sur Particle 
de 1933 ne fait pas non plus mention de Aitken. 

2. Le compte-rendu de Zentralblatt est similaire. Dans les deux cas, il n’est pas question du procede d’Aitkeii. 
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paraitre surprenant, mais nous aliens voir que les conditions de vie et de travail de Georges 
Lemaitre en cette annee de guerre 1942 ainsi que ses preoccupations scientifiques expliquent 
cela. 

Mots cles Lemaitre, Aitken, Steffensen, A^, point fixe, equation differentielle ordinaire, 
methode iterative de Picard. 
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1 Le de Aitken et la methode de Steffensen 


1.1 L’article de Aitken 


Le resume de Particle de AitkenH est le suivant : 


Proceedings 

XXV-On Bernoulli’s Numerical Solution of Algebraic Equations 
A. C. Aitken 

The aim of the present paper is to extend Daniel Bernoulli’s method of 
approximating to the numerically greatest root of an algebraic equation. On the 
basis of the extension here given it now becomes possible to make Bernoulli’s 
method a means of evaluating not merely the greatest root, but all the roots of an 
equation, whether real, complex, or repeated, by an arithmetical process well adapted 
to mechanical computation, and without any preliminary determination of the nature 
or position of the roots. In particular, the evaluation of complex roots is extremely 
simple, whatever the number of pairs of such roots. There is also a way of deriving 
from a sequence of approximations to a root successive sequences of ever-increasing 
rapidity of convergence. 


Ainsi le point de depart du travail de Aitken est le memoire de Daniel Bernoulli (1700-1782) 
Observationes de seriebus quae formantur ex additione vel subtractione quacunque terminorum 
se mutuo consequentium \BeTno\illi D., 1728| . A une equation algebrique, mise sous une forme 
particuliere, Bernoulli fait correspondre une suite recurrente lineaire. En utilisant une ecriture 
moderne, si Pequation est 1 = aix -\- ■ —h OpX^, la suite associee est Un+p = aiUn+p-i H—h ctpUn 
> 0)0. Avec de bonnes hypotheses, le rapport —— tend vers la solution de plus petit 


[n 


^n+1 


module de Pequation. Une analyse de cette methode est donnee dans [Chabert J.-L. and, 1994 


pages 254-259. Les auteurs indiquent que Leonhard Euler dans [Euler L., 1748 donne une 
preuve de la methode de Bernoulli en considerant une fraction rationnelle de denominateur 
1 — aix — ... — OpX^ qu’il developpe en serie entiere de deux fagons. Tout d’abord, il fait 
remarquer que les coefficients du developpement suivent la relation de recurrence0. Puis il 
obtient le developpement en decomposant la fraction rationnelle en elements simples en faisant 
done apparaitre les racines de Pequation. Euler d’ecrit d’ailleurs le cas d’une equation a racines 


simples : Pelement simple 


1 


1 — pz 


se developpe sous la foYmeJ2k>o ^ ce qui permet d’avoir 


pages 


la regie de Bernoulli0. On retrouve toute cette demarche exposee dans [Henrici P., 1988 
584-588. 

Aitken0 commence par generaliser la methode precedente de Daniel Bernoulli dite methode 
des series recurrentes. Partant d’une equation -p bp_iz^~^ -f • • • -P biz -P 5o = 0 et de la suite 


3. Le resume est disponible sur : http ://dx.doi.org/10.1017/S0370164600022070 

4. C’est d’ailleurs la presentation choisie par Henrici [Henrici R, 1988| 

5. a partir d’un certain rang. 

6. Euler indique de plus que pour obteiiir la racine de plus grand module de P{z), on applique la methode 

au polynome P etant de degre p. 

7. Une analyse de son article est faite dans [Chabert J.-L. and, 1994} pages 492-497. 


4 





















recurrente Sn+p + bp_iSn-\-p-i + • • • + + b^Sn = 0, il considere les rapports : 


'^n+l 


^n-\-2 

•^n+l 


'^n+m 
'^n+m—1 


'^n—m+2 

^n—1 


S 


n—m+3 


'^n+l 


^n+1 


'^n+m—1 
'^n+m—2 


( 1 ) 


'^n—m+1 '^n—m+2 


'^n 


Aitken prouve que la limite de ce rapport quand n ^ +oc est Z 1 Z 2 - - - oil Zi, Z 2 ^ • • •, Zm sont 
les m premieres racines de Tequation, supposees telles que > \z 2 \ > • • • > 12:^1. 

En permutant les lignes dans chacun des determinants, on obtient le rapport : 


'^71—m+2 '^n—m+3 


'^n+1 




'^n+l 


'^n+m—1 


•^n+l 

m+1 


^n+2 


'^n+m 


'^n—m+2 


( 2 ) 


'^n—1 

^n+1 


'^n+m—2 
'^n+m—1 


qui est en fait le quotient de deux determinants de Hankel (1839-1873) associes a la serie 

j^n—m+2 

J2n>o^nZ^ ■ ■ Rappelons que I’expression generale des determinants de Hankel est : 




"Si+l 

^i+i — 1 

H] = 

•^z+l 

Si+2 • 

Si+j 



Si+j 

• Si+2j-2 


Peter Henrici donne une demonstration de ce rfeultat de convergence dans [Henrici P., 1988| , 
pages 596-6021^. II lie ensuite ce resultat a I’algorithme quotient-difference de Heinz Rutishauser 
(1918-1970) IRutishauser H., 1954| . 


1.2 Acceleration de la convergence 

Aitken cherche ensuite a ameliorer la vitesse de convergence de suites qu’il a prealablement 
considerees. Ce qu’il fait en considerant les differences finies des premier et second ordre de la 

suite {tn = ]. En utilisant son resultat ci-dessus pour m = 2, il montre que les differences 

V Sn / 

8. Il fait d’ailleurs reference a Particle de Aitken. 
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premieres se comportent comme les termes d’une suite geometrique de raison —. Aitken pose 
alors la relation ; 

in+2 _ ^^n+1 


^1 ^ n +1 

En resolvant par rapport a Zi, il vientl^ : 


Atr. 


Zl = 


tntn+2 tn+1 

AH„. 


(5) 


Le terme de droite est ce que I’on nomme aujourd’liui le transforme par le d’Aitken de la 
suite (tn). On espere que cette nouvelle suite va converger vers zi plus vite que la suite initiale. 
Aitken propose ensuite de faire subir a cette nouvelle suite la transformation du pour avoir 
une convergence plus rapide. Comme le rappelle Claude Brezinski dans [Brezinski C., 1996 


Aitken indique que la transformation du A^ a d^a ete consideree par d’autres auteurs mais 
pas dans un but d’acceleration de la convergence. On pent nuancer un pen cela, toujours 
dans [Brezinski C., 199^ (voir aussi [Brezinski C., 19^ et [Osada N., 2012[ ), Claude Brezinski 
indique que le mathematicien japonais Takakazu Seki (7-1708) introduit la transformation 
du A2 vers 1680 dans le calcul du volume d’une sphere. D’apres Osada, la methode de Seki 
correspond a une acceleration de convergence!^. 


1.3 La methode de Steffensen 


L’article de Steffensen [Steffensen, J. F., 1933[ est publie en 1933 dans les Skandinavisk 
Aktuarietidskrift, une revue destinee aux actuaireso dont il a ete un des editeurs 
pendant trente ans [Ogborn M.E., 1962 . Contrairement a Aitken et a Lemaitre, le point 

un memoire original d’un mathematicien celebre. 
The calculus of observations de Whittaker- 


de depart de Steffensen n’est pas 
Steffensen fait reference a deux ouvrages 
Robinson [Whittaker E.T. and Robinson G., 1924 et Vorlesungen uher numerisches Rechnen 
de Runge-Konig. Steffensen considere une equation scalaire f{x) = x dont il cherche a approcher 
une racine reelle. Il debute son article par quelques considerations sur la suite des iterees 
Xn+i = f{xn) selon si / est croissante on decroissante (sur un intervalle). Ceci lui permet 
d’encadrer la racine. Steffensen introduit la notion de suites iterees opposees. 

Ensuite Steffensen raisonne en termes d’interpolation polynomiale. Il interpole la fonction 
/ aux points (iterfe) xq, Xi = f{xi) et X 2 = /(a;i)@en utilisant les differences diviseesEl : 

X = f{x) = f{xo) + (x - Xo) [f{,Xo),f{Xi)] + (x - Xo)(x - xi) [/(x), /(xo), /(xi)] (6) 


9. Dans Particle de Aitken, le nnmeratenr de la fraction qni suit est ecrit sous la forme d’un determinant 
2x2. 

10. Osada ecrit : « From the numerical analysis view point, the formula is the m-panels trapezoidal rule for 

2 

the numerical integration / Ay{D — y)dy = Therefore we can say that Seki aceelerated the trapezoidal 

Jo 3 

rule by the Aitken ISJ process. » 

11. Steffensen a debute sa carriere dans les assurances a Copenhague. 

12. Comme I’indique Steffensen, tout choix de trois points consecutifs convient pour I’etude. 

13. voir par exemple [Davis P.J., 1975| pages 39-40 et 64-65. 
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ou encore 


X2 — Xi 


(7) 


X ^ Xi + {x — Xq) 


Xi — Xq 


Cette derniere expression condnit a : 



( 8 ) 


Steffensen developpe alors plusieurs exemples nnmeriqnes reposant sur I’algorithme dont nous 
donnons une ecriture classique : 



Sur les exemples nnmeriqnes proposes par Steffensen, I’amelioration de la vitesse de convergence 
est bien visible. 


Si deux exemples proposes par Steffensen portent sur des calculs d’interets, Steffensen devait 


etre aussi motive par des questions de Dynamique de populations. En effet, en 1930, Steffensen 


reprend une question sur le probleme d’extinctioneiudiee alors par Agner Krarup Erlang (1878- 


1929) [Bacaer N., 2011 . On revient toujours a une recherche de point hxe. 


2 L’iteration rationnelle de Georges Lemaitre 


2.1 L’annee 1942 


Quel est le contexte historique et quelles sont les preoccupations scientihques de Georges 


Lemaitre en 1942 ? Nous nous appuyons ici essentiellement sur la biographie de Dominique 


Lambert [Lambert D., 200'0| . En Mai 1940, la Belgique est envahie par I’armee allemande. La 
bibliotheque de I’Universite Catholique de Louvain est detruite par le feu comme en Aout 
19140 Apres avoir tente de gagner I’Angleterre avec plusieurs membres de sa famille, Georges 
Lemaitre revient en Belgique en Juin 1940. II demenage alors une serie de machines a calculs 
mecaniques de type . Georges Lemaitre utilise ces machines pour ses calculs et 

vent eviter de plus qu’elles ne tombent entre les mains de I’occupant. Georges Lemaitre est 
sensibilise an moment de ses etudes universitaires a cette question du calcul sur machine par 



ses professeurs (puis collegues), Edouard Goedseels (1857-1928) et surtout Maurice Alliaume 



|. A partir des annees 40, Georges Lemaitre va consacrer une part importante 


de son activite scientihque an calcul sur machines. Le point d’orgue est I’introduction -et la 
programmation- du premier ordinateur a I’Universite Gatholique de Louvain en 19580. Dans 


Jones G.B., Lloyd J.L. (eds)201l|, Pierre-Jacques Gourtois dans sa contribution The Belgian 


Electronic Mathematical Machine (1951-1962) : An Account, indique que Gharles Manneback 

14. site de I’Universite Catholique de Louvain 

15. Lambert op.cit., page 218. 

16. Lambert op.cit., pages 42-43. 

17. voir aussi la biographie de Georges Lemaitre sur le site de I’UCL. 
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(1894-1975) a porte le projet de construction d’un ordinateur beige an debut des annees 50. 
Or Manneback, professeur de Physique-Matheniatique et d’Electricite a I’ecole d’ingenieurs de 
Louvain, est un collegue et ami de Georges Lemaitre. Dans les annees cinquante, en Suisse, un 
projet du meme type voit le jour. Le mathematicien Heinz Rutishauser, cite plus haut et qui a 
generalise -et implemente- la metliode de Aitken, participe a la creation du premier ordinateur 
suisse, ERMETH [Bauer F.L., 2002 . 


2.2 La methode 


Georges Lemaitre propose une modification de la methode des approximations successives, 
Xn+i = f{xn) pour rfeoudre I’equation x = /(x)@ Notons que Lemaitre fait reference 
comme Steffensen aux deux ouvrages : The calculus of observations de Whittaker- 
Robinson [Whittaker E.T. and Robinson G., 1924 


et Vorlesungen iiber numerisches Rechnen 
de Runge-Konig. Deux ouvrages classiques de I’epoque, qui devaient se trouver etre facilement 
accessibles. A I’instar de Aitken, Georges Lemaitre construit son propos a partir des travaux 
d’un mathematicien illustre, en I’occurence Garl Friedrich Gauss (1777-1855). Geci n’est pas 
surprenant. Dominique Lambert!^ souligne que la lecture dans les textes originaux (en latin) 
faisait partie des habitudes de Lemaitre, habitude heritee de Henri Bosmans (1852-1928)o 
D’ailleurs Lambert precise que Georges Lemaitre prive de ses contacts internationaux durant la 
Seconde Guerre mondiale, s’est plonge en particulier dans la Mecanique Analytique de Lagrange 
(1736-1813) et dans Les methodes nouvelles de la mecanique celeste de Henri Poincare (1854- 
1912). 

Georges Lemaitre adapte a I’iteration x = f{x) la mGhode de Gauss. Gauss s’interesse a 
un systeme 

j X{x,y) = 0 . . 

in^,i/) = o 


dont il connait trois valeurs approchees {xn,yn), n = 1,2,3. II s’agit de dGerminer a partir 
des ces trois valeurs, une valeur pour {x,y). Gauss linearise le probleme en considerant les 
developpements de Taylor a I’ordre 1, sans reste, an point {x,y) des functions X et Y. En 
Waluant ces deux developpements en {xn, Vn) pour n = 1,2, 3, on obtient trois systemes lineaires 
d’inconnues respectives (x — x„,?/ — yn) {n = 1,2,3). Gauss donne ensuite I’expression pour 
un n fixe, de x — Xn en fonction de Xn = X{xn,yn), de W = Y{xn,yn) et des coefficients des 
developpements de Taylor. Geux-ci etant inconnus. Gauss ecrit une condition de compatibilite 
entre les trois expressions sous la forme d’un determinant nul. Georges Lemaitre reprend cette 


18. Lemaitre indique que cette methode a ete etendue a la resolution d’equations differentielles par Emile 
Picard (1856-1941). C’est justement la motivation de Lemaitre comme on le voit d’ailleurs a la fin de son 
article. 

19. Lambert op.cit, page 222. 

20. Concernant la vie et I’oeuvre du pere jesuite Henri Bosmans, on pent se referer a 
[Hermans M., Stoffel J.-F., (eds.)2010] . 
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demarche. II linearise /, 


{ ^2 = f{xi) = a + /3xi 
X3 = f{x2) = a + 13x2 
X = f{x) = a + j3x 

puis donne la condition de compatibilite du systeme : 


X2 

1 

Xi 

X3 

1 

X2 

X 

1 

X 


( 11 ) 


( 12 ) 


En developpant le determinant selon la derniere ligne, on obtient : 

X1X3 — Xn 

X = - 

Xi - 2 X 2 + X3 


(13) 


c’est a dire la transformation du d’Aitken! 

Si Georges Lemaitre ne donne pas d’exemple numerique, il considere bien la methode comme 
une methode d’acceleration de la convergence. En effet, il conclut : 

On pent done employer la formule d’iteration rationnelle pour accelerer la 
convergence la ou elle existe. 


Cette phrase suit les considerations faites par Georges Lemaitre sur la resolution d’une equation 

differentielleE^ = f{x,t) pour laquelle on connait une solution approchee X 2 {t). En suivant 
dt 

le processus iteratif de Picard [Picard E., 1890| , pages 197-200 


X3 


f{x 2 ,t)dt ; 


(14) 


Lemaitre propose alors d’appliquer la methode d’iteration rationnelle aux trois fonctions xi{t), 
X 2 {t) et X 2 ,{t). L’objectif de Georges Lemaitre est en fait la question de la resolution approchee 
d’equations differentielles. Il publie dans le meme numero du Bulletin de la Classe Scientifique 
un article sur I’application de la methode d’iteration rationnelle a I’integration d’une equation 
differentielle ordinaire [Lemaitre G., 1942 b| dont nous allons parler. 


2.3 Un second article de Georges Lemaitre 


Goerges Lemaitre propose done dans [Lemaitre G., 1942 b 


une methode d’integration 
numerique pour la resolution d’equations differentielles ordinaires. Il developpe son propos 
en traitant I’exemple de I’equation 


s' = ^ = - x) 


(15) 


on y est fonction de x. Il s’agit pour lui de calculer des valeurs approchees de la solution y 
qui tend vers x quand x tend vers -|-cx). Georges Lemaitre n’explique pas pourquoi une telle 


21. probleme avec condition initiale. 
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solution existe. Georges Lemaitre va cependant considerer I’isocline 1 (1 = — x)) et la 

droite y = x (I’isocline 0). En fait, il y a un anti-entonnoir compris entre I’isocline 1 et la droite. 
Le signe de la derivee partielle de y'^{y — x) par ra pp ort a y permet de conclure qu’il y a une 
solution et une seule piegee dans cet anti-entonnoirri. 

Nous avons essaye de voir ce qui se passe avec un ordinateur. Maplel^, le logiciel choisi, utilise 
une metliode de Runge-Kutta-Fehlberg (RKF45)0. Nous n’avons pas pu mettre en evidence 
cette solution en travaillant dans le sens positif, i.e. sur I’intervalle [0,+cxo[. Cependant, on 
trouve des valeurs semblables jusqu’a 2.2 a cedes trouvees par Lemaitre et rfeumees dans la 
tableau se trouvant a la fin de Particle (il calcule les valeurs sur I’intervalle [0; 3]. La valeur de la 
condition initiale donnee par Lemaitre, x(0) = 0.618343 joue bien un role de « separateur »dans 
le sens oil si I’on prend une valeur legerement superieure, la solution sort de I’anti-entonnoir par 
le haut, tandis que si I’on prend une valeur inferieure, la solution traverse I’isocline 0. On pent 
mettre en evidence la solution en travaillant dans le sens negatif, i.e. sur I’intervalle ]—cxo, 0] car 
la solution est alors stable dans ce sens0 

Georges Lemaitre n’explique pas pourquoi il a choisi cette equation^. On pent imaginer 
qu’il ait pu choisir cette equation parce que les methodes classiques d’integration numerique 
ne fonctionnent pas de fagon satisfaisante. Il applique done la methode decrite a la fin de 
[Lemaitre G., 1942 a : il part d’une courbe liee a I’equation differentielle qui donne la « premiere 
approximation yi », il derive et integre yi pour obtenir deux nouvelles functions yo et y 2 - Georges 
Lemaitre applique ensuite I’iteration rationnelle aux valeurs yoix), yi{x) et y 2 {x). Pour choisir 
la function initiale yi, Lemaitre considere la famille des isoclines p = 2y‘^{y — x) dont il forme 
I’equation differentielle associee : 


, ^ ^ ^ y 

^ dx 3y — 2x' 


(16) 


Par des conditions geometriques (qualitatives), Lemaitre considere la « courbe limite »donnee 
par la condition ; 



y 

3y — 2x 


p = 


(17) 


soit 


2yi(vi -x)(3pi -2x) = l 


(18) 


courbe somme toute proche de I’isocline 1. Il lui reste a la parametrer. 

22. Pour toutes ces notions ont peu se referer a [Hubbard J., West B., 199^ . 

23. sur la version 14 de Maple 

24. qui est une methode explicite. 

25. Nous remercions le professeur Ernst Hairer de nous avoir signale cela et d’avoir fait des essais numeriques 
en utilisant son code RADAU5 (Runge-Kutta implicite). En partant d’une valeur initiale y{x*) = x*, il obtient 
dans la direction negative les valeurs suivante pour j/(0) : si a;* = 1.0, 1/(0) = 0.602680285037; si x* = 3.0, 
1/(0) = 0.618340077402 ; si x* = lOh i/(0) = 0.618340077404; si x* = lOh y{0) = 0.618340077404; si x* = 10^ , 
1/(0) = 0.618340077404; si x* = 10®, i/(0) = 0.618340077404. La precision est de 3 * 10“^® et le nombre de pas 
d’integration est de 582 pour x* = 10®. 

26. Georges Lemaitre s’interesse particulierement a cette epoque au probleme de Stormer et developpe des 
algorithmes pour les machines mecaniques de calcul numerique qu’il utilise a Louvain. 
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La methode de Georges Lemaitre est dont un melange de methodes qualitatives (choix de 
yi), de la methode iterative de Picard (avec aussi des methodes pour les calculs annexes) et du 
d’Aitken pour ameliorer les resultats. 


3 Conclusion 


L’objectif de cette note n’etait pas de regler une question de primaute d’un resultat 
scientifique. Les exemples abondent en mathematiques de methodes redecouvertes puis 
amGiorees ou de rfeultats prouvfe independamment. Les preoccupations des auteurs et les 
dynamiques de recherches ont leur importance dans ce phenomene. Dans le cas de I’iteration 
rationnelle, c’est la resolution de certaines equations differentielles qui incite Georges Lemaitre 
a chercher un moyen d’accelerer une methode connue (celle de Picard) et qui soit facile a mettre 
en oeuvre, en un mot a implementer sur une machine a calculs mecanique. 

Les motivations des trois mathematiciens etaient bien distinctes. Aitken se demarque, son 
objectif de depart est le calcul approche des racines d’une equation polynomiale a une variable. 
Steffensen et Lemaitre se rejoignent : il s’agit d’accelerer la convergence d’un processus iteratif. 
Mais Steffensen en reste au stade des iterees d’une function a une variable, tandis que Lemaitre 
applique sa methode a des functions de plusieurs variables, mais surtout a des equations 
differentielles 

L’isolement scientifique et les affres de la guerre expliquent en partie que Lemaitre n’ait pas 
eu connaissance de I’existence des travaux de Aitken et de Steffensen. Ajoutons que Lemaitre 
avait une fagon bien a lui de travailler, en citant Lambert^ : 


Remarquons aussi un detail qui est caracteristique de la maniere dont Lemaitre 
travaillait. II poursuit les idees qui I’interessent en ne se souciant guere de la 
litterature existante ou des groupes qui pourraient ahorder les memes sujets que 
lui. 


27. Lambert op.cit., page 147 
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4 Annexes 


4.1 Compte-rendu du Jahrbuch de Particle de Aitken 

Aitken, A. C. 

On Bernoulli’s numerical solution of algebraic equations, Proceedings Royal Soc. 
Edinburgh 46, 289-305. 

Die Bernoullische Methode zur numerischen Auflosung einer algebraischen 
Gleichung 

g{x) = + ■ ■ • + = 0 

beruht auf folgendem Satze : Wird bei beliebigen Anfangsgliedern durch die 
Rekursionsgleichung 


aofit + n) + aif{t + n - 1) -\ -h Onfit) = 0 


die Folge f{t) bestimmt, so ist 


lim 

t—>-OD 


.fit + l) 

fit) 


= zi, 


wobei Zi die Wurzel vom grofiten absoluten Betrage bedeutet; allerdings mufi noch 
vorausgesetzt werden, dafi die andern Wurzeln sdmtlich von kleinerem Betrage sind 
als Zi. 

Diese Bernoullische Formel verallgemeinert nun Verf. dahingehend, dafi die 
Mdglichkeit gegeben wird, sdmtliche Wurzeln gleichzeitig zu bestimmen. 

Setzt man ndmlich fm{t) gleich, der Determinante 


\f{t + X-K)\ 


so ist 


lim 

t^OO 


fmjt + 1 ) 

fmit) 


= 2 ^ 12^2 


wobei Zi, Z 2 , ■ ■ ■, Zn die Wurzeln von g{x) = 0 in absteigender Ordnung nach ihrem 
absoluten Betrage bedeuten und vorausgesetzt wird, dafi\zm\ > |^m+i| sei. 

Die Praxis bietet auch fiir die Falle, in denen die allgemeine Formel versagt, 
Mdglichkeiten, die Wurzeln einer algebraischen Gleichung nach der Regel zu 
bestimmen, die durch den angegebenen Satz angedeutet ist. Diese Mdglichkeiten 
werden in den weiteren Abschnitten der Arbeit an Hand einiger Beispiele diskutiert. 

Im weiteren Verlauf der Untersuchungen werden aus dem Vorangehenden noch 
Methoden hergeleitet, schdrfere Approximationen der Wurzeln zu erzielen. 

Zum Schlufi wird die Methode des Verf. noch mit andern verglichen, insbesondere 
mit der Methode von Fiirstenau (Darstellung der reellen Wurzeln algebraischer 
Gleichungen durch Determinanten der Goefficienten, Marburg 1860), die als ein 
Spezialfall (bei passender Wahl der Anfangselemente der Folge fit)) erscheint. (IV 
17.) 

W.; Dr. (Konigsberg i. Pr.) 
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4.2 Compte-rendu du Jahrbuch de Particle de Steffensen 

Steffensen, J. F. 

Remark on iteration, Skand. Aktuarietidskrift 16, 64-72,1933 

Wird die zu losende Gleichung auf die Form x = f{x) gehraucht, und setzt man 
= f{xy), so ergibt sich ein guter Anndherungswert aus einer Folge dreier Werte 
Xo,Xi,X 2 durch 

(Axo)^ 

^ = ^0 - 7 ^-• 

Lorey, W.; Prof. (Frankfurt am Main) 


4.3 Compte-rendu du Jahrbuch de Particle de Lemaitre sur Piteration 
rationnelle 

Anschliefiend an eine von Gaufi in der Theoria motus entwickelte Methode 
gibt Verf. ein Verfahren an, aus mehreren durch einfache Iteration gewonnenen 
Ndherungswerten einen neuen besseren Wert zu gewinnen. Fiir eine Gleichung mit 
einer Unbekannten ergibt sich unter der Annahme, dafi die Funktion in dem in Frage 
kommenden Gebiet linear angendhert werden darf, aus der Vertrdglichkeitsbedingung 
der drei Gleichungen 

^2 = fi^i) = a + bxi, X 3 = f{x 2 ) = a + bx 2 , x = f{x) = a + hx 
der neue Ndherungswert 

X1X3 — xl 

X = -. 

Xi - 2X2 + X 3 

Der Vorteil dieser Formel ist der, dafi man auch zum Ziel kommt, wenn die Folge 
der Ndherungswerte xi, X 2 , x^ nicht gegen x konvergiert (das Verfahren stimmt mit 
einem von Steffensen, Skand. Aktuarietidskr. 16 (1933), 64-72 (F. d. M. 58, 535) 
gegebenen uberein). Bei Gleichungen mit mehreren Unbekannten erhdlt man fiir jede 
der Unbekannten eine entsprechende Gleichung. 

Zum Schlufi wird kurz auf die Anwendung des Verfahrens auf Differentialglei- 
chungen eingegangen. Dabei unterscheidet Verf. eine Iteration durch Integration und 
eine durch Differentiation. 

Willers, F.; Prof. (Dresden) 


4.4 Compte-rendu du Jahrbuch de Particle de Lemaitre sur la 
resolution d’une equation differentielle par Piteration rationnelle 

Integration d’une equation differentielle par iteration rationnelle. 

Acad. Belgique, Bull. Cl. Sci. (5) 28, 815-825. 

Das in der vorstehend besprochenen Arbeit beschriebene Verfahren wird auf die Gleichung 

0 2/ N 

dx 
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angewendet. Aus einer ersten Naherung yi werden zwei weitere durch 


— = 2 / (yi - x)dx und yo - x = - - 

y2 ^ dx yi 

gewonnen und aus diesen in der angegebenen Art eine neue Naherung berechnet. Gesucht wird 
eine Losung im Intervall 0 < x < cxd, die die Gerade y = x als Asymptote hat. Ausgegangen 
wird von der Wendepunktskurve der Differentialkurven 

‘2yi{yi - x){3yi - 2x) = 1 , 

die zunachst punktweise mit dem Argumentsschritt Ax = 0,1 von x = 0 bis x = 3 berechnet 
wird. Dazu wird yo gebildet und berechnet und weiter durch numerische Integration y 2 . Die aus 
diesen drei Naherungen berechnete neue ist auf sechs Dezimalen genau. Ftir grohe x wird nach 
der Methode der unbestimmten Koeffizienten eine asymptotische Reihe bestimmt, aus der sich 
ftir X = 3 bis auf sechs Stellen der gleiche Wert ergibt wie der oben berechnete. Angabe der 
Resultate in Tabellen. 

Willers, F.; Prof. (Dresden) 
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